Depinning transition in failure of inhomogeneous brittle materials 
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The dynamics of a crack propagating in an elastic inhomogeneous material is investigated. The 
variations of the average crack velocity with the external loading are measured for a brittle rock 

and are shown to display two distinct regimes: Below a given threshold G c , the crack velocity is 

c 

well described by an exponential law v ~ e G ~( r > characteristic of subcritical propagation, while 
for larger values of the driving force G > G c , the velocity evolves as a power law v ~ (G — G c )° 
with 8 — 0.80 ± 0.15. These results can be explained extending the continuum theory of Fracture 
Mechanics to disordered systems. In this description, the motion of a crack is analogue to the one of 
an elastic line driven in a random medium and critical failure occurs when the loading is sufficiently 
large to depinne the crack front from the heterogeneities of the material. 

PACS numbers: 62.20.Mk, 46.50.+a, 68.35. Ct 



Failure of inhomogeneous materials has been a very ac- 
tive field of research during the last decades (see Ref. [l| 
for a recent review). A great research effort in this field 
has been dedicated to the study of fluctuations: Fluc- 
tuations of velocity around the average motion of cracks 
when the studies were devoted to their highly intermit- 
tent dynamics @, S II H, or variations to a straight 
trajectory when the works were dedicated to the rough 
geometry of fracture surfaces (EBB- In both cases, 
these fluctuations were shown to display remarkably ro- 
bust properties suggesting that crack propagation in dis- 
ordered systems could be described on a general manner 
by relatively simple statistical models able to capture the 
competition between the two antagonist effects occurring 
during failure of inhomogeneous materials: Disorder and 
elasticity. 

Very recently, main statistical features of fluctuations 
of both trajectory and velocity for cracks propagating in 
brittle materials were captured by stochastic models of 
elastic lines driven in random media 0, [Toj ] that mimic 
the motion of cracks through the microstructural disorder 
of materials. However, the relevance of this theoretical 
framework for fracture problems is still a matter of de- 
bate: On the one hand, the ability of these models to de- 
scribe the average behavior of the crack such as its mean 
velocity, or the critical external loading at failure, more 
interesting from a mechanical or an engineering point of 
view, is still an open question. On the other hand, a 
direct experimental observation of the critical dynamic 
transition from a crack pinned by the heterogeneities of 
the material (v = 0) to a propagating crack (v > 0), as 
predicted by this theory at the onset of material failure 
(driving force G = G c ), is still lacking. The investigation 
of this depinning transition on an experimental example 
is the central point of this Letter. 

The variations of the average crack velocity with the 
external driving force, i.e. the energy release rate G 111 ]. 



are measured for a brittle rock. They are shown to ex- 
hibit two distinct regimes. Below a critical threshold G c , 

the crack velocity is well described by an exponential law 

_ c 

v ~ e G -< r > characteristic of a subcritical propagation, 
while for larger values of the external loading G > G c , 
the velocity evolves as a power law v ~ (G — G c ) with 
6 = 0.80 ± 0.15. This behavior is fully captured by 
a stochastic model directly derived from Fracture Me- 
chanics extended to inhomogeneous systems where crack 
propagation is analogue to the motion of an elastic line 
driven in a random medium 1^, T3, 13j [- 
System and setup. - 

Sandstone is chosen as an archetype of heterogeneous 
elastic materials. A Botucatu sandstone, extracted in 
the central region of Brazil, has been used for the ex- 
periments. It is made of quartz grains with a diameter 
d = 230 /jm ± 30 /xm and a porosity <f> = 18 ± 2 %, 
that results in highly inhomogeneous mechanical prop- 
erties at the grain scale. This South-American rock is 
consolidated thanks to an iron oxide cement providing 
to the rock a red coloration. As a result, its fracture 
energy G c ~ 140 J.m~ 2 as measured in the following is 
relatively high compared to other sandstones Its 
intrinsic tensile strength measured by splitting cylinders 
submitted to uniaxial compression [15J is found to be 
ay = 75 MPa ± 5 MPa, leading to an estimate of its pro- 
cess zone size £pz — § ^r- [16( - zone next to the crack 
tip where damage and dissipative processes are localized 
during failure - of £pz — 200 /im. This failure mecha- 
nisms property (£pz < d) suggests that crack propaga- 
tion in the Botucatu rock mimics the crack motion in 
ideal brittle materials at the scale of the quartz grains. 

A new experimental setup has been developed in order 
to measure variations of crack velocity from slow to very 
fast propagation in brittle materials. Contrary to vari- 
ous fracture tests as the Double Torsion geometry [T3| 
classically used to measure loading/crack velocity rela- 



FIG. 1: Experimental setup, (a) Sketch of the Tapered Dou- 
ble Cantilever Beam geometry; (b) picture of the specimen 
during crack propagation. 
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tionship in rocks, the Tapered Double Cantilever Beam 
specimens used in the experiments (Fig.QJa)) result in 
a slight, but controlled acceleration of the crack that is 
produced by the tapered shape of the specimen. As a 
consequence, it is possible to measure crack velocities up 
to v ~ 1 ms' 1 not achieved by classical tests used in 
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|. In addition, failure is obtained 



by the propagation of a straight crack front in the spec- 
imen, so that the local velocity at each position of the 
crack line is also the average velocity of the front, allow- 
ing for great simplifications in the theoretical analysis of 
the experiments. Finally, a straight crack propagation 
in the specimen is obtained without any lateral guide 
grooves that are known to induce systematic errors on 
experimental v(G) curves (l8j . 

An initial notch cq — 35 mm is machined in 100 mm 
long samples with thickness e = 30 mm. They are then 
submitted to a uniaxial traction by increasing the dis- 
placement <5f = Vcxtt between two rods inserted in the 
drilled specimens at constant velocity 0.2 /im s _1 < 
v C xt < 4 /im s" 1 . During the test, a gauge force mea- 
sures the applied tension F while a clip gage measures 
the opening displacement 5 between the two lips of the 
crack with a precision of 0.1 /im (see Fig.QJb)). A typical 
force-crack opening displacement curve obtained during 
a fracture test of a Botucatu specimen at room temper- 
ature is presented in Fig. [21 

The initial linear part of the curve - prior crack initi- 
ation - allows for an estimation of the Young's modulus 
E = 25 ± 1 GPa of the sandstone, in agreement with 
the value obtained from measurements of its compressive 
and shear waves speed. After crack initiation, the average 
position of the crack front c = (c(z)) z is measured using 
Finite Element (FE) simulations of an elastic specimen in 
the same geometry. The numerical compliance A FE (c) is 
compared to the experimental compliance X(t) = S/F in 
order to get the evolution of the crack length c(t) repre- 
sented in inset of Fig.[2j Other techniques limited to the 
free surface of the sample based either on image analysis 
of the crack motion at the surface or on the resistance 
measurement of a thin conductive film deposited on the 



FIG. 2: Mechanical behavior of the specimen, (a) Typi- 
cal load-crack opening displacement curve; (b) corresponding 
evolution of the average position of the crack front. 



sample side have led to similar, however less precise mea- 
surements of the crack length. 

From the evolution of the crack length, it is now pos- 
sible to measure the crack speed v = 4| as well as the 
driving force G imposed to the system during the test. 
Using the load-displacement curve to measure the work 
SW of the tensile machine during the span <5i[l||, one gets 

G W = elc(ilTi)- C (t)j B- 0n the 0ther hand ' the driv " 

ing force is estimated independently using the relation 
G(t) = [F(t))] 2 g FE [c(t)} where the geometrical part g FE 
of the driving force is provided by the FE simulations. 
Both methods lead to similar results within 2 %. 

Experimental results. - The variations of the crack ve- 
locity with the driving force as observed on the sandstone 
specimens are represented in Fig. [3] in semi-logarithmic 
coordinates. Velocity measurements are achieved over 
almost five orders of magnitude, corresponding to a rel- 
atively small variation of the driving force. Irrespective 
of the external loading rate as well as the sample in- 
vestigated, the failure behavior of the rock is found to 
be systematically characterized by two distinct regimes 
defining G c . Near, but above this critical loading, a slight 
change in the driving force results in a strong variation 
in the crack velocity. This high sensibility is studied in 
more detail in the bottom right inset of Fig.[3j where 
v is plotted as a function of the distance to the criti- 
cal force G — G c in logarithmic coordinates. The linear 
behavior in this representation suggests a power law vari- 
ation of the crack velocity v ~ (G — G c ) e . The value of 
G c = 140 ±3 J.m~ 2 is found to optimize this scaling rela- 
tion, and leads to an exponent 9 — 0.80 ±0.15 where the 
error bars are calculated from the variations measured 
from sample to sample. 

The variations of velocity at low driving forces G < 
G c are now studied. Contrary to the previous regime, 
slow crack propagation in rocks has been largely inves- 
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tigated and shown to depend crucially on the tempera- 
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Analytical forms as v 

[22| are usually used to describe the experi- 
mental data. Both formula reproduce correctly the mea- 
surements reported here as far as G < 120 J.m~ 2 . The 
first one, used by many researchers because it is conve- 
nient to integrate and differentiate leads to a subcritical 
crack growth index n ~ 34 that compares well with the 
other experimental findings for sandstone 14]. The sec- 
ond formula leads to b ~ 0.68 . 10~ 20 m 2 which is also in 
agreement with the other measurements made on rocks 
with a similar structure [l8j]. This description is based 

on the Arrhenius law v ~ e fc B T where the activation en- 
ergy E a — E — bG represents the typical barrier along 
the energy landscape of the system tilted by the external 
force G. In this approach, the tensile force produced by 
the external loading on the interatomic bond just next 
to the crack tip favours the thermal activated processes 
responsible for its rupture, as e.g. thermal stress fluctua- 
tions [2^] and chemical reactions [24| . As a consequence, 
the motion of the crack front is analogue to the prop- 
agation of one point in a 2D medium. As soon as one 
considers a 3D inhomogeneous material, the elasticity of 
the crack line comes into play, and one expects a different 
expression of the activation energy E a ~ g-{t) 1^' ^ as 
presented in detail in the next section. This expression 
provides a rather good description of the data for the full 
subcritical regime G < G c . The best linear regression of 
the function ln(u) ~ y>(— gjrf) ) represented in the upper 
left inset of Fig.[3]is obtained for (r) = 63 ± 5 J.m- 2 . 

Discussion. - The observation of two very different 
regimes, with an exponential variation of v with the ex- 
ternal loading for G < G c and a power law behavior 
for G > G c , reveals a fundamental aspect of the dy- 
namics of crack propagating in brittle inhomogeneous 
medium. Let's investigate theoretically the crack mo- 
tion in these systems to understand quantitatively such 
a behavior. As a starting point, we assume that the local 
velocity v(M) of a point M of the crack is proportional 
to the excess of energy G(M) — T(M) locally released by 
the system where T refers to the local fracture energy 
of the material. In a disordered material such as sand- 
stone, the fracture energy can be described as a stochas- 
tic field r(M) = (r) + STr](M) where 77 is a short range 
correlated random term with zero mean value and unit 
second order moment. The heterogeneities of the ma- 
terial induce perturbations of the crack front, which is 
parallel in average to the z-axis and propagates along 
the x-axis, both in the mean fracture plane (x, z) (in- 
plane perturbations c(z, t) — (c(z, t)) z ) and in the perpen- 
dicular direction y (out-of-plane perturbations h(z,t)). 
They in turn lead to variations in the local value of 
the external driving force G(M). Interestingly, for small 
perturbations, G(M) is only depending on the in-plane 




FIG. 3: Average dynamics of a crack propagating in Botucatu 
sandstone. The variations of the crack velocity are plotted in 
logarithmic scale with respect to the external loading. The 
subcritical regime G < G c with G c = 140 J.m -2 is studied in 
the top-left inset. Solid line corresponds to the best fit of the 

data in v ~ e G < r > obtained for (r) = 63 J.m -2 . Bottom 
right inset shows the velocity variations with the net loading 
G — G c in a logarithm representation for G > G c . Straight 
line corresponds to a power law fit with exponent 8 = 0.81. 



deviations of the crack front [26j, 

G f c(z,t)-c(z',t) ^ 



, and is given by 
where G refers to 



the macroscopic driving force applied by the tensile ma- 
chine to the specimen. Using the previous expressions of 
the local driving force and material fracture energy, one 
gets the following equation of motion for a crack propa- 
gating in a brittle inhomogeneous material 



§l,.~«Mr»4 



c(z,t)-c(z',t) 
(z-z') 2 



dz +8Tr)(c, h, z). 

(1) 

As the out-of-plane perturbations h behave indepen- 
dently of c, the stochastic term in this equation is ana- 
logue to a 2D random potential depending only on c 
and z. Therefore, the crack motion is described by an 
equation of pinning/dcpinning of an elastic line driven 
in a random medium as previously proposed for inter- 
facial cracks propagating in inhomogeneous weak planes 
[l2l Il3j |: if the external loading G, or driving force for 
unit crack length, exceeds a given threshold 



G c = (r) + ^T) 2 /(r), 



(2) 



the crack propagates, while the front is pinned by the ma- 
terial heterogeneities if G < G c [29]. Above the thresh- 
old, the mean velocity v of the crack front is expected to 
scale as (G — G c ) where 9 is called the velocity expo- 
nent. Ertas and Kardar have studied equation ([T]) using 
functional renormalization group technique [3(3] ■ To first 
order in perturbation, they find 9 = 0.78. Recent direct 
numerical simulations resulted in 9 ~ 0.63 1311 . As a con- 
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sequence, the power law behavior measured experimen- 
tally with an exponent 9 ~ 0.80 suggests that a depinning 
transition from a stable to a propagating crack as de- 
scribed in Eq. |T]) occurs at G — G c . Below the threshold 
at zero temperature, the external driving force is not suf- 
ficient to make the crack propagate and the crack front is 
pinned by the material heterogeneities. However, at finite 
temperature T > 0, thermally activated processes can en- 
able a subcritical propagation. In this regime, described 
by adding an annealed noise rf^{z,t) to Eq. ([I]), one ex- 
pects also a collective motion of the line, characteristic 
of glassy systems, and velocity variations are predicted 

fi) 2 <r> 2 

to follow the so-called creep law v ~ e. fe B T < G -< r >) 



33l | . This expression describes rather well the ex- 



perimental measurements of v(G) over the whole range 
of external loading G < G c , leading to activation ener- 
gies in the range E a = (j) 2 ~ 0.20 - 0.30 eV, 
one order of magnitude larger than the thermal energy 
k^T . Interestingly, the expression of the activation en- 
ergy provides an estimate of the topothesy - size of 
the elementary feature of the crack front for which the 
amplitude of its roughness compares with its length - 
£ = 0.10 nm, compatible with the interatomic distance. 
Finally, let us note that using Eq. @, the experimental 
values of the critical loading G c and the average fracture 
energy (r) allows to estimate the normalized fluctuations 

0.62 ± 0.06 of fracture energy in the 



sr 



-<r> 



^<r> 



Botucatu rock, slightly larger but comparable with an 

0.5 corresponding to 



estimate of this quantity y 
an ideal porous material made of an homogeneous solid 
with constant fracture energy and voids, with volumic 
fractions 1 — <j> = 0.82 and <fi = 0.18, respectively. 

Conclusion. - The average dynamics of a crack propa- 
gating in a brittle inhomogeneous rock has been investi- 
gated. The velocity variations with the external loading 
display two distinct regimes: above a threshold G c , v 
evolves as a power law of the net loading G — G c with 
exponent 8 = 0.80 ± 0.15 while for G < G c , these varia- 
tions are described by an Arrhenius law v ~ e fc s T with 
an activation energy E a ~ gz^r) ■ Both behaviors can be 
explained extending the continuum theory of fracture me- 
chanics to inhomogeneous systems. In this description, 
the external loading, the crack front elasticity as well as 
the effect of material heterogeneities interplay leading to 
a macroscopic motion of the line above a critical force 
G c larger than the average value of the fracture energy 
in the material (r). Below this critical transition, the line 
can propagate only at finite temperature through ther- 
mal activated processes. Therefore, the measurement at 
the macroscopic scale of the crack velocity variations with 
the external loading represents a signature of the depin- 
ning of the crack front at the scale of the material mi- 
crostructure. The experimental measurements reported 
here represent a strong argument in favour of the pin- 



ning/depinning approach to brittle failure and open new 
perspectives for the prediction of macroscopic quantities 
such as average crack velocity and fracture energy of di- 
rect interest for Engineering and Applied Science. 
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